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Abstract 

The  Role  of  Ceiling  Points  in 
General  Integer  Linear  Programming 

\ 

\  Robert  M.  Saltzman  and  Frederick  S.  Hillier 

\  Stanford  University,  1988 

This  report  examines  the  role  played  by  several  kinds  of  ceiling  points  in  solving  the  pure, 
general  integer  linear  programming  problem  (ILP).  While  no  assumptions  are  made  concerning 
the  structure  or  signs  of  the  data  of  the  problem,  it  is  assumed  that  the  feasible  region  for  (ILP) 
is  non-empty  and  bounded.  A  ceiling  point  with  respect  to  a  single  constraint  may  be  thought  of 
as  an  integer  solution  on  or  close  to  the  boundary  of  the  feasible  region  defined  by  the  constraint. 
The  definition  of  a  ceiling  point  with  respect  to  a  single  constraint  is  extended  to  take  multiple 
constraints  into  consideration  simultaneously,  defining  what  is  called  a  feasible  ceiling  point.  It  is 
shown  that  the  set  of  all  feasible  ceiling  points  contains  at  least  one  optimal  solution  for  (ILP).  A 
related  class  of  solutions  called  feasible  1 -ceiling  points  is  also  characterized  and  shown  to  contain 
all  optimal  solutions  for  (ILP).  Moreover,  1 -ceiling  points  are  computationally  easier  to  identify 
than  ordinary  ceiling  points  and  may  be  sought  with  respect  to  one  constaint  at  a  time.  It  is  also 
demonstrated  that  solving  (ILP)  requires  only  enumerating  feasible  1-ceiling  points  with  respect 

y  * 

to  a  subset  of  all  functional  constraints.  ■  -  k-  ,•  ■->  *'  > 
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1.  Introduction  and  Overview 


In  this  report  we  consider  the  role  played  by  ceiling  points  in  solving  the  pure,  general 
integer  linear  programming  problem  ( ILP )  in  m  constraints  and  n  variables  x y,  j  = 
whose  form  is 

•  •  T 

maximize  c  x  =  z 

subject  to  Ax  <  b  (ILP) 

x  >  0,  x  integer, 

where  A  e  9?"**",  b  £  9?m  and  c  6  9in.  All  the  data  {A,6,c}  are  assumed  to  be  rational 
numbers,  but  they  are  unrestricted  in  sign.  The  problem  is  “pure”  in  that  all  of  the 
variables  are  required  to  take  on  nonnegative  integer  values.  It  is  “general”  in  the  sense 
that  the  variables  may  take  on  any  nonnegative  integer  values  permitted  by  Ax  <  b,  as 
opposed  to  being  restricted  to  0  or  1  (the  binary  case).  An  important  additional  assumption 
is  that  no  implicit  or  explicit  equality  constraints  are  used  to  define  the  feasible  region 
FR  =  {x  >  0|  Ax  <  b).  This  implies  that  any  m'  x  n'  subsystem  of  equality  constraints, 
where  m'  <  n',  has  already  been  reexpressed  by  solving  for  m'  of  the  variables  in  terms  of 
the  other  n'  —  m!  variables  and  substituting  for  these  m'  variables  in  all  of  the  remaining 
inequality  constraints,  yielding  a  problem  consisting  only  of  inequality  constraints. 

Our  goal  is  to  develop  a  solid  understanding  of  the  nature  and  utility  of  ceiling  points. 
In  the  next  section,  ceiling  points  are  defined  and  some  of  their  properties  explored.  Section 
3  extends  these  basic  definitions  to  those  of  d-ceiling  points,  which  are  computationally 
easier  to  identify  than  ordinary  ceiling  points.  In  Section  4,  we  consider  what  constraints 
are  most  critical  in  the  search  for  an  optimal  solution,  while  the  rcstrictiveness  of  assump¬ 
tions  made  in  various  theorems  is  examined  in  Section  5.  Section  6  summarizes  the  key 
concepts  discussed  in  this  report  using  an  example  in  the  plane. 
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2.  Characterizing  Ceiling  Points 


Before  formalizing  the  intuitive  concept  of  a  ceiling  point  as  an  integer  solution  on  or 
close  to  the  surface  of  the  feasible  region  for  the  linear  programming  relaxation  ( LPr )  of 
(ILP),  a  few  conventions  are  noted.  Throughout  the  report,  constraint  (i)  may  refer 
to  a  functional  or  nonnegativity  constraint  of  the  form  ajy  <  If  it  refers  to  the  jth 
nonnegativity  constraint,  all  coefficients  a,j  are  0,  except  for  the  jth  which  is  —1,  and  the 
corresponding  right  hand  side  bi  is  0.  Also,  the  term  “direction”  will  be  used  to  mean  a 
direction  in  3R"  which  is  parallel  to  one  of  the  n  coordinate  axes.  The  first  ceiling  point 
definition  considers  an  integer  point  satisfying  a  single  constraint  i  with  little  or  no  slack. 

Definition  1  An  integer  vector  a;  is  a  ceiling  point  with  respect  to  the  ith  constraint, 
denoted  x  =  CP(i),  if 

(1)  of x  <  bi,  and 

(2)  ajx  +  la^l  >  bi,  for  each  j  =  1,  ...,n. 

The  first  part  of  the  definition  simply  means  that  x  satisfies  the  ith  constraint;  the 
second  part  of  the  definition,  called  the  ceiling  point  condition,  implies  that  taking  a  unit 
step  from  x  toward  the  ith  constraining  hyperplane  in  every  direction  results  in  an  infeasible 
point.  Letting  xj  =  x  +  e7  and  xj  =  x  —  e;  represent  the  points  which  are  a  unit  step 
away  from  x  in  the  plus  j  and  minus  j  directions,  respectively,  we  can  reexpress  (2)  from 
our  definition  as  follows.  If  >  0  then  xj  violates  constraint  (*),  whereas  if  ai;  <  0  then 
xj  violates  constraint  (i),  for  all  j.  If  atJ  =  0  then  changing  the  jth  component  of  x  does 
not  affect  the  feasibility  of  x  with  respect  to  constraint  (t).  f 

|  It  may  be  interesting  to  compare  Definition  1  with  that  first  given  by  J.  Wilson  [Wi, 
Definition  1],  Assuming  that  atl  >  al2  >  . . .  >  ajn  >  0,  6,  >  0,  and  all  data  axe  integer, 
he  defines  a  nonnegative  integer  vector  x  to  be  a  ceiling  point  with  respect  to  a  single 
constraint  afy  <bi  if 

(1)  aTx  < 
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Thus,  according  to  the  ceiling  point  condition  (2)  of  the  definition,  no  ceiling  points 
with  respect  to  (t)  alone  exist  when  any  =  0.  This  suggests  the  need  for  a  second  defi¬ 
nition  which  considers  not  just  one  constraint,  but  all  of  the  nonnegativity  and  functional 
constraints  of  ( ILP )  simultaneously. 

Definition  2  An  integer  vector  a:  is  a  ceiling  point  with  respect  to  the  feasible  region 
FR  =  {x|  Ax  <  b,x  >  0},  denoted  x  =  CP(FR),  if 

(1)  x  €  FR,  and 

(2)  3 *  :  afx  +  |ajj)  >  for  each  j  =  1,  ...,n. 

In  this  case,  x  =  CP(FR)  indicates  that  x  is  feasible  and  for  every  direction  j  either 
x+  or  xj  (or  both)  violate  at  least  one  functional  or  nonnegativity  constraint.  Notice  that 
x  =  CP(i )  and  x  €  FR  imply  that  x  =  CP(FR).  However,  the  converse  is  not  true, 
i.e.,  x  =  CP(FR)  does  not  imply  that  x  =  CP(i),  for  some  (i),  as  illustrated  by  point  x 
in  Figure  1. 

Ceiling  points  in  an  integer  linear  program  play  a  role  analogous  to  that  of  corner-point 
feasible  solutions  in  a  linear  program. 


Lemma  la  Every  extreme  point  of  the  convex  hull  of  feasible  integer  solutions  for  {ILP)  is 
a  CP(FR ). 


Proof:  First  note  that  all  extreme  points  of  the  convex  hull  of  feasible  integer  solutions 
are  themselves  feasible  integer  solutions  (since  the  vertices  of  the  convex  hull  of  a  finite  set 


(2)  afx  +  ay  >  bi  +  1,  for  j  =  1, ...,  n,  and 

(3)  Xj+1  >  0  =>  afx  +  a,  -  aJ+i  >  6j  +  1,  1  <  j  <  n. 

In  contrast,  Definition  1  makes  no  assumptions  on  the  relative  magnitudes,  signs  or  inte¬ 
grality  of  the  coefficients,  nor  does  it  require  the  vector  x  to  be  nonnegative.  Wilson’s  goal 
is  to  first  locate  n  ceiling  points  satisfying  the  ith  constraint  with  equality,  and  then  use 
these  points  to  construct  tne  minimal  inequality  corresponding  to  this  constraint,  i.e.,  to 
strengthen  this  linear  integer  inequality  by  reducing  its  coefficients  as  much  as  possible. 
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U  *  CP(FR); 


V  =  CP(  1); 


W  =  CPC  2); 

X  =  CP(  FR); 

V  =  CP(  2); 


u  »  CP(  i),  i  =  1 ,2.3 


V  =  CP(FP) 

w*  CP(FR) 


x  *  CP(  i),  i  =  1 ,2,3 


y  =  CP(FR) 


Figure  1.  Examples  of  Ceiling  Points. 


of  points  are  a  subset  of  the  points  from  which  the  convex  hull  is  formed  [PS,  p.  37]).  Now 
suppose  that  x  is  such  an  extreme  point  which  is  not  a  CP(FR).  Then  there  is  a  direction 
j  such  that  both  xj  and  xj  are  feasible.  Hence,  x  =  (x+  -f  xj  )/2,  contradicting  the 
extreme  nature  of  x.  I 


Because  they  tend  to  be  closer  to  the  boundary  of  FR  than  non-ceiling  points,  ceiling 
points  in  the  “interesting”  portion  of  FR  are  likely  candidates  for  having  the  highest 
objective  function  values.  This  idea  has  motivated  our  interest  in  ceiling  points:  they 
seem  to  be  the  points  on  which  we  should  focus  our  attention.  Indeed,  we  can  show  that 
if  an  optimal  solution  x*  for  ( ILP )  exists,  either  x*  itself  is  a  CP(FR)  or  there  exists 
another  optimal  solution  which  is  a  CP(FR). 


2.  Characterizing  Ceiling  Points _ 5 

Theorem  1  Suppose  the  set  of  feasible  solutions  for  ( ILP )  is  non-empty  and  bounded. 
Then  there  exists  an  optimal  solution  for  (ILP)  which  is  a  CP(FR). 

Proof:  Since  (ILP)  possesses  at  least  one  feasible  integer  solution  and  is  bounded,  it 
possesses  an  optimal  solution.  Furthermore,  the  convex  hull  of  feasible  integer  solutions  is 
non-empty.  By  a  well-known  theorem  (see  [HL,  p.  54]  or  [Si,  p.  91])  an  optimal  solution  for 
a  linear  programming  problem  occurs  at  one  of  the  extreme  points  of  its  feasible  region. 
Similarly,  in  an  integer  linear  program,  an  optimal  solution  occurs  at  an  extreme  point  of 
the  convex  hull  of  feasible  integer  solutions.  Since  all  such  extreme  points  are  CP(FR)' s, 
by  Lemma  la,  it  follows  that  at  least  one  optimal  solution  for  (ILP)  is  a  feasible  ceiling 
point,  a  CP(FR).  I 

The  analogy  between  ceiling  points  and  extreme  points,  however,  should  not  be  car¬ 
ried  too  far.  First  notice  that  not  every  ceiling  point  is  an  extreme  point  of  the  convex 
hull  of  feasible  integer  solutions  for  ( ILP)  as  evidenced  by  the  ceiling  point  v  in  Figure  1 . 
Moreover,  while  only  a  finite  number  of  extreme  points  in  the  feasible  region  of  ( LPr )  may 
exist,  the  feasible  region  of  the  associated  (ILP)  may  possess  an  infinite  number  of  ceiling 
points  if  it  is  unbounded.  Even  if  the  feasible  region  is  bounded,  the  number  of  CP(FR)' s 
is  likely  to  be  much  larger  than  the  number  of  extreme  points  in  the  associated  (LPr). 
Nonetheless,  the  main  implication  of  the  theorem  is  that  an  (ILP)  may  be  solved  by  enu¬ 
merating  and  comparing  the  values  of  all  CP(FR)' s:  the  ceiling  point(s)  with  the  largest 
objective  value  optimally  solve(s)  the  problem.  Also,  if  a  complete  enumeration  is  not 
possible,  it  might  be  useful  to  concentrate  one’s  effort  on  locating  ceiling  points  relatively 
near  x,  where  any  feasible  solution  found  is  apt  to  be  a  high-quality  and  sometimes  even 
optimal  solution.  Locating  such  ceiling  points  is  the  goal  of  the  Heuristic  Ceiling  Point 
Algorithm  described  in  a  subsequent  report. 

Corollary  lb  Suppose  the  set  of  feasible  solutions  to  (ILP)  is  non-empty  and  bounded. 
Further  suppose  all  components  of  c  differ  from  zero.  Then  every  optimal  solution  for 
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( ILP )  is  a  CP(FR). 

Proof:  Since  (ILP)  possesses  at  least  one  feasible  integer  solution,  it  also  possesses  an 
optimal  solution.  Let  x*  be  an  optimal  solution  for  (ILP)  which  is  not  a  CP(FR).  Then, 
since  x*  is  feasible,  there  must  exist  at  least  one  direction  j  such  that  both  x*+  and  x*~ 
are  feasible.  Consider  any  such  direction  j.  Since  Cj  ^  0, Vj,  either  x*+  or  x*~  has  a 
strictly  superior  objective  value  to  that  of  x*,  (because  either  cTx*+  =  cTx*  +  cTe.}  = 
cTx*  +  cj  >  cTx *,  or  cTx*~  =a  crx*  —  cTc.j  =  cTx *  —  Cj  >  crx*),  contradicting  the 
optimality  of  x*.  I 

A  few  more  observations  can  be  made  about  the  utility  of  ceiling  points,  or  lack  thereof. 
When  some  c}  —  0,  the  proof  of  Corollary  lb  implies  that  there  may  exist  non-ceiling 
points  which  are  optimal.  On  the  other  hand,  if  the  optimal  solution  for  (ILP)  is  unique, 
it  must  be  a  CP(FR)  by  Theorem  1.  However,  we  cannot  assume  a  priori  that  the  optimal 
solution  is  unique  nor  is  it  realistic  to  assume  that  all  c}  differ  from  0.  For  instance,  in  fixed- 
charge  problems,  the  fixed-charge  (0-1)  variables  sometimes  have  a  zero  cost  component. 
Consequently,  we  have  found  it  beneficial  to  relax  the  ceiling  point  definitions  given  above. 


3.  d-Ceiling  Points  and  Their  Properties 


From  a  computational  standpoint,  identifying  a  CP(FR)  is  difficult:  in  order  to 
satisfy  the  ceiling  point  condition,  (part  (2)  of  the  definition),  a  constraint  must  be  violated 
as  a  result  of  either  increasing  or  decreasing  every  component  by  one  unit.  We  now 
describe  two  sets  of  points  which  are  formed  by  relaxing  the  ceiling  point  conditions  of 
Definitions  1  and  2,  respectively,  so  that  a  constraint  violation  does  not  have  to  occur  in 
every  direction  from  the  ceiling  point.  Consequently,  these  sets  are  easier  to  enumerate 
than  their  respective  ceiling  point  counterparts. 
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Definition  3  For  an  integer  d  €  [1,  n],  an  integer  vector  x  is  a  d-ceiling  point  with  respect 
to  the  ith  constraint,  denoted  x  =  d-CP(i),  if 

(1)  afx<  and 

(2)  afx  +  \a,ij\  >  bj ,  for  each  of  d  distinct  indices  j  —  ji , jj. 

In  other  words,  x  =  d-CP(i)  means  x  satisfies  the  ith  constraint,  and  taking  a  unit 
step  from  x  toward  the  ith  constraining  hyperplane  in  directions  j  =  j\ ,  ...,jd  (and  perhaps 
in  some  additional  directions)  results  in  an  infeasible  point.  For  example,  if  atjl  >  0  then 

violates  constraint  (i),  whereas  if  a,jt  <  0  then  x~  violates  constraint  (t).  As  with 
ordinary  ceiling  points,  Definition  3  may  be  extended  to  cover  a  set  of  constraints  which 
comprise  the  feasible  region  FR  of  an  ( ILP ). 

Definition  4  For  an  integer  d  €  [l,n],  an  integer  vector  x  is  a  d-ceiling  point  with  respect 
to  the  feasible  region  FR,  denoted  x  —  d-CP(FR),  if 

(1)  x  €  FR,  and 

(2)  3 i  :  ajx  +  |o,;|  >  b,,  for  each  of  d  distinct  indices  j  = 

In  this  case,  x  =  d-CP(FR)  indicates  that  x  is  feasible  and  in  directions  j  =  j i, ...,  jj 
either  x~j  or  xj  ( or  both)  violate  at  least  one  functional  or  nonnegativity  constraint.  It 
should  be  clear  that  a  point  which  is  a  d-CP{i)  for  a  particular  d  is  also  a  d-CP(i)  for 
every  integer  d  €  [l,d  —  1],  The  same  is  true  for  d-CP(FR)' s.  The  points  in  Figure  1 
may  be  recast  as  d-CP(i)' s  as  follows:  v  =  2-CP(l),  w  =  2-CP(2),  x  =  l-CP(l)  and 
x  =  1-CP(2),  while  y  =  2-CP(2)  and  y  =  1-CP(3).  It  is  worth  emphasizing  that  a  single 
point  may  be  a  d-CP(i)  with  respect  to  more  than  one  constraint,  especially  if  d  is  allowed 
to  vary.  As  with  ceiling  points,  x  =  d-CP(i)  and  x  €  FR  imply  that  x  =  d-CP(FR). 
Again  referring  to  Figure  1,  the  points  v,x,  and  y  all  are  d-CP(FR)' s,  with  d  =  2. 

We  now  wish  to  clarify  the  relationships  between  the  various  sets  of  ceiling  points 
introduced.  The  definitions  indicate  that  an  n-CP(i)  is  equivalent  to  a  CP(i ),  as  is  an 


3.  d-  Ceiling  Points  and  Their  Properties 


8 


n-CP(FR)  to  a  CP(FR).  However,  when  d  <  n,  the  ceiling  point  condition  for  a  d- 
CP(i)  is  less  restrictive  than  that  for  a  CP(i),  as  it  is  for  a  d-CP(FR)  in  comparison  to 
that  for  a  CP(FR).  This  implies  that  the  cardinality  of  the  set  of  ordinary  ceiling  points 
with  respect  to  a  constraint  ( i )  is  no  larger  than  the  cardinality  of  the  set  of  d-ceiling 
points  with  respect  to  that  same  constraint,  for  any  specific  integer  d  £  [l,n].  Just  as  the 
set  of  CP(i)'s  with  respect  to  a  particular  (i)  is  a  subset  of  the  set  of  that  constraint’s 
d-CP(i)’ s,  the  set  of  CP(FRYs  is  a  subset  of  the  set  of  all  d-CP(FRys.  The  following 
lemma  summarizes  some  of  the  key  relationships  between  the  four  types  of  ceiling  points, 
and  is  illustrated  for  d  =  1  by  Figure  2. 

Lemma  lc  For  any  integer  d  £  [l,n]  the  following  relationships  exist  among  the  4  sets  of 
ceiling  points  identified  by  Definitions  1  -  4: 

(a)  {*|  x  =  CP(FR)}  C  {r|  x  =  d-CP(FR)} 

(b)  U,  {*1  *  =  CP(i)}  C  U.  {*|  x  =  d-CP(i)} 

(c)  {x|  r  =  d-CP(FR)}  C  U,  {x|  x  =  l-CP(i)} 

(d)  (J,  {*1  *  =  feasible  d-CP(i)}  C  {x|  x  =  d-CP(FR)} 

T  oof:  (a)  From  the  ceiling  point  conditions  of  Definitions  2  and  4,  we  have  {x|  x  =  1- 
CP(FR)}  D  {x\  x  =  2 -CP(FR)}  D  .  .  .D  {x|  x  =  n-CP(FR)}  =  {x|  x  =  CP(FR)} 
Thus,  every  CP(FR)  is  a  d-CP{FR).  Figure  2  depicts  this  for  the  case  where  d  =  1:  the 
set  of  CP(FRY s,  labelled  as  box  3,  is  competelv  contained  within  the  set  of  l-CP(FJ?)’s, 
labelled  as  box  4. 

(b)  From  the  ceiling  point  conditions  of  Definitions  1  and  3,  {x|  x  =  1-CP(?)}  D 
{x|  x  =  2 -CP(i)}  2  •••  5  {x|  x  —  n-CP(i)}  =  {x|  x  =  CP(i)}.  Since  this  holds  foi  all 
constraints  individually,  it  also  holds  for  the  union  over  all  constraints.  This  is  shown  for 
d  =  1  in  Figure  2  by  the  box  labelled  2  being  completely  contained  within  box  1,  the  union 
of  the  set  of  l-CP(i)’s. 
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(c)  Let  x  be  any  d-CP(FR).  Then  there  is  at  least  one  direction  such  that  either 
x+  or  xj  violates  some  constraint.  Thus,  x  is  a  l-CP(i)  for  at  least  one  (i).  However, 
equality  in  (c)  does  not  hold  since  a  l-CP(i)  may  violate  a  constraint  other  than  (t).  With 
d—  1,  the  box  labelled  4  in  Figure  2  is  completely  contained  within  box  1. 

(d)  Let  x  be  any  feasible  d-CP(i).  Then  x  satisfies  the  definition  of  d-CP(FR). 
However,  equality  in  (d)  does  not  necessarily  hold  for  d  >  1  because  a  CP(FR)  need  not 
be  a  d-CP(i)  with  respect  to  any  particular  constraint  (as  illustrated  for  d  —  2  by  the 
point  x  in  Figure  1).  With  d  =  n,  the  union  over  (i)  of  the  set  of  all  feasible  n-CP(iys  is 
shown  in  Figure  2  as  the  intersection  of  boxes  2  and  3,  labelled  box  5.  I 


Feasible  Region 


{x|  x=CP(FR)} 


(x|  X  =  l-CP(FR)} 


U  {x|  x  =  CP ( 1)} 

i 


Y 


{x|  x  =  1  -CP( 0} 


{x|  all  integer  points  in  tne  unit  byperspnere  about  x  are  feasible} 


Figure  2.  Relationships  Between  the  4  Types  of  Ceiling  Points  (Lemma  lc). 


In  general,  the  larger  the  magnitude  of  d  in  d-CP(i),  the  closer  a  d-ceiling  point 
is  to  the  surface  of  the  feasible  region.  However,  an  n-CP(i)  is  not  necessarily  a  more 
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desirable  solution  than  a  1  -CP(i).  First,  the  n-CP(i)  may  be  infeasible  while  the  1- 
CP(i)  is  feasible.  Second,  assuming  feasibility,  the  objective  function  value  of  these  two 
integer  solutions  depends  more  upon  their  location  in  the  feasible  region  relative  to  the 
objective  function  hyperplane  than  it  does  upon  d.  What  is  important  to  note  is  that  even 
a  feasible  l-CP(i),  i.e.,  a  1  -CP(FR),  may  be  an  optimal  solution,  as  demonstrated  in  the 
next  theorem. 

Theorem  2  Suppose  the  set  of  feasible  solutions  for  ( ILP )  is  non-empty  and  bounded. 
Let  Sc  =  {j|  c}  0}  be  the  support  of  c.  Then  every  optimal  solution  for  (ILP)  is  a 
d-CP(FR),  where  d  >  |SC|. 

Proof:  Since  (ILP)  possesses  at  least  one  feasible  integer  solution  and  is  bounded,  it 
also  possesses  an  optimal  solution.  Let  x*  be  an  optimal  solution  which  is  not  a  |SC|- 
CP(FR).  Then,  for  at  least  one  j  €  Sc,  both  x*+  and  x*~  are  feasible.  However,  since 
Cj  ^  0 ,Vj  e  Sc,  either  x*+  or  x*~  has  a  strictly  superior  objective  value  to  that  of  x*, 
contradicting  the  optimality  of  x*.  I 

In  particular,  if  |5C|  =  1,  then  all  optimal  solutions  for  (ILP)  are  l-CP(FJ?)’s.  This 
fact  should  not  be  too  surprising  since  {xj  x  =  l-CP(FF)}  excludes  only  those  points 
which  are  “interior”  to  the  feasible  region  in  the  sense  that  all  integer  points  within  the  unit 
hypersphere  about  them  are  feasible.  (These  are  the  points  in  the  box  labelled  0  in  Figure 
2  which  are  not  in  any  other  box.)  While  the  set  of  CP(FR)' s  contains  the  set  of  vertices 
of  the  convex  hull  of  all  feasible  integer  solutions  to  (ILP),  the  set  of  \-CP(FR)'s  contains 
all  integer  points  on  the  surface  of  this  convex  hull,  including  its  vertices.  Consequently, 
another  way  of  solving  (ILP)  is  to  enumerate  the  feasible  elements  of  (J(  {x|  x  =  1- 
CP(i)}  because  among  them  are  all  the  l-CP(PP)’s  for  (ILP).  The  advantage  of  taking 
this  route  is  that  we  can  focus  on  finding  l-CP(i)’s  with  respect  to  one  constraint  at  a 
time,  while  checking  simultaneously  for  feasibility.  This  is  the  central  approach  of  the 
algorithms  presented  in  subsequent  reports  for  solving  an  integer  linear  program. 
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When  |SC|  >  1,  can  we  solve  ( ILP )  by  enumerating  all  feasible  |Sc|-CP(i)’s  with 
respect  to  one  constraint  at  a  time?  Unfortunately,  if  jSc|  >  1  and  we  search  only  for 
feasible  we  may  fail  to  find  some  \Sc\-CP(FR)'s,  by  Lemma  lc(d).  Moreover, 

it  is  possible  that  the  only  optimal  solution  for  an  (ILP)  is  a  feasible  l-CP(t)  which  is 
not  also  a  feasible  d-ceiling  point  (d  >  1)  with  respect  to  any  single  constraint.  In  this 
case,  enumerating  only  feasible  |Sc|-CP(t)’s  will  fail  to  yield  the  optimal  solution.  For 
example,  consider  the  following  simple  (ILP): 

max{xi  +  x2  +  X3I Xj  <  4.5,.;  =  l,2,3;x  >  0}. 

The  integer  solution  y  =  (4,4,4)  is  a  3 -CP(FR),  yet  y  is  a  l-CP(i)  with  respect  to  each 
upper  bound  constraint  individually.  In  fact,  no  2-CP(v)’s  or  3-CP(i)’s  exist.  Here,  y  is 
the  unique  optimal  solution,  and  it  would  not  be  found  by  searching  only  for  |5'c|-CP(?)’s. 
Thus,  if  we  take  the  approach  of  seeking  d-CP(i)' s  with  respect  to  one  constraint  at  a 
time,  we  must  take  d  =  1. 


4.  Search  Constraints 


Suppose  we  have  an  algorithm  G  for  finding  l-CP(t)’s  with  respect  to  a  given  con¬ 
straint  (i).  Then  Theorem  2  indicates  that  to  solve  (ILP)  we  need  to  apply  our  algorithm 
G  to  each  functional  and  nonnegativity  constraint.  The  purpose  of  this  section  is  to 
demonstrate  that  we  only  need  to  search  for  l-CP(t)’s  with  respect  to  a  subset  of  all  the 
constraints.  Also,  this  set  of  “search  constraints”  is  not  difficult  to  identify.  The  next 
theorem  indicates  that  we  may  exclude  from  the  set  of  search  constraints  all  of  the  non¬ 
negativity  constraints  because  if  there  exists  an  optimal  solution  which  is  a  1-ceiling  point 
with  respect  to  a  nonnegativity  constraint,  there  also  exists  an  optimal  solution  which  is 
a  1-ceiling  point  with  respect  to  a  functional  constraint. 


4.  Search  Constraints 


12 


Theorem  3  Suppose  the  set  of  feasible  solutions  for  (ILP)  is  non-empty  and  bounded. 
Further  suppose  that  (ILP)  is  not  trivial  to  solve  (by  trivial  we  mean  that  if  x  =  0  is 
feasible  then  it  is  optimal).  Then  (ILP)  possesses  an  optimal  solution  which  is  a  feasible 
1- CP(i),  where  (i)  is  a  functional  constraint. 

Proof:  Let  x*  be  any  optimal  solution  for  (ILP).  Assume  that  neither  x*  nor  any  other 
optimal  solution  is  a  l-CP(i),  where  (t)  is  some  functional  constraint. 

For  all  j  =  l,...,n  there  are  three  possibilities: 

Case  I.  Cj  >  0  :  Neighboring  solution  x*+  violates  no  functional  constraints  since  x*  ^  1- 
CP(FR),  nor  does  it  violate  any  nonnegativity  constraints  because  x*  >  0.  Furthermore, 
crx*+  >  cTx* ,  contradicting  the  optimality  of  x*. 

Case  II.  cj  =  0  :  Solutions  along  a  ray  emanating  from  x*  all  have  the  same  objective 
value,  i.e.,  cTx*  =  cT(x*  +  tj)  =  ...  =  cT(x*  -f  Kej).  Eventually,  for  some  I<  large 
enough,  x*  +  I\ej  violates  a  functional  constraint  (fy)  by  the  boundedness  assumption. 
Then  x*  +  (K  —  l)ey  is  a  feasible,  optimal  1  -CP(ij),  contradicting  the  hypothesis. 

Case  III.  Cj  <  0  :  If  neighboring  solution  x*~  is  feasible,  then  cTx*~  >  cTx*,  contra¬ 
dicting  the  optimality  of  x*.  Thus,  the  only  remaining  possibility  is  that  x*~  is  infeasible, 
implying  that  Xy  =  0,  which  is  considered  next. 

Since  these  cases  hold  for  every  direction  j,  we  are  left  with  x*  =  0  for  all  j,  i.e.,  x*  =  0 
is  the  only  optimal  solution  for  an  (ILP)  with  c  <  0.  But  this  is  precisely  the  situ¬ 
ation  ruled  out  by  the  theorem’s  suppositions:  If  x  =  0  is  a  feasible  solution  for  the 
(ILP)  max{cTx|  Ax  <  b,  x  >  0,  c  <  0},  then  x  =  0  must  be  optimal.  Therefore,  ei¬ 
ther  x*  itself  or  another  optimal  solution  is  a  l-CP(i)  with  respect  to  some  functional 
constraint  (i).  I 

We  now  demonstrate  that,  at  least  in  R2,  we  need  not  search  for  l-CP(i)'s  with 
respect  to  any  constraint  that  does  not  intersect  the  unit  cube  with  all-integer  vertices 
containing  x,  an  optimal  solution  for  (LPr).  To  help  show  this,  let  A  =  {t|  afx  =  6,}  be 
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the  set  of  constraints  binding  at  x,  and  C\  =  {x|  o(x  <  6, ,  Vi  £  A)  be  the  cone  formed 
by  the  extreme  rays  of  FR  emanating  from  x.  Also  let  K  =  {j |  Xj  =  integer}  be  the  set 
of  indices  of  integer-valued  components  of  x,  k  =  |A'|  and  UHC[x]  be  any  of  the  2“  unit 
hypercubes  in  with  all-integer  vertices  which  contain  x.  If  k  =  0,  i.e.,  xj  ^  integer, 
Vj,  then  there  is  a  unique  UHC\x\  =  {x6  [x^J  <  xj  <  Jxj"],V;}.  If  k  =  1,  i.e.,  x/  = 
integer,  for  some  /,  then  there  are  two  corresponding  UHC[x ]’s:  {x  £  3in|  xj  —  1  <  x/  < 
xt;  [xj\  <  xj  <  [xj],Vj  ^  i}  emd  {x  £  Sfin|  xt  <  xt  <  xj  +  1;  |xjj  <  xj  <  fx;] , V;  /  /}. 
If  k  1,  then  there  exist  2*  UHC[x\ s  defined  in  the  obvious  manner.  Since  it  is  assumed 
that  x  is  not  all-integer,  the  scalar  k  £  {0,1,..., n  —  1}. 

Theorem  4  Suppose  an  ( ILP )  with  n  =  2  is  non-empty  and  bounded.  Further  suppose 
that  c  ^  0  and  x  is  an  optimal  solution  for  (LPr).  Then  every  optimal  solution  for 
{ILP)  is  a  feasible  l-CP(i),  where  constraint  (i)  intersects  each  of  the  2*  UHC[x]'s. 

Proof:  The  proof  will  be  by  contradiction  in  a  rather  lengthy  case  analysis.  In  each 
case,  we  first  assume  the  existence  of  an  optimal  solution  y  G  for  {ILP)  which  is  not 
a  1  -CP{i)  with  respect  to  any  (i)  intersecting  a  UHC[x],  We  then  reach  a  contradiction 
by  specifying  a  nonnegative  integer  solution  x*  which  is  both  feasible  and  strictly  better 
than  y  in  any  possible  configuration  of  x,  the  binding  constraints,  y  and  c  satisfying  the 
original  assumptions.  Note  that  if  c  =  0  then  all  feasible  integer  solutions  are  optimal,  not 
just  1-CP(FjR)’s,  so  the  theorem  does  not  hold. 

We  first  define  a  region  which,  under  the  assumptions  about  y,  must  contain  only 
feasible  integer  solutions.  The  general  picture  to  keep  in  mind  is  shown  in  Figure  3(a).  The 
points  {ui,  U2,  «3,  U4}  are  the  all-integer  vertices  of  UHC[x],  while  the  points  {iq,  v2,  v3,  vt } 
are  the  all-integer  vertices  of  what  Balas,  et  a  1.  [BB]  would  refer  to  as  the  “dual  to  the  unit 
hypercube  centered  at  y,”  abbreviated  here  as  DU HC[y\.  Figure  3(a)  shows  a  “feasibility 
cone”  Ci  whose  vertex  is  at  x  and  whose  edges  are  formed  by  constraints  binding  at  x.  The 
edges  of  Ci  that  will  be  used  in  the  proof  and  that  are  shown  in  Figure  3(a)  are  the  tightest 
possible  “admissable  binding  constraints,”  i.e.,  constraints  binding  at  x  which  satisfy  the 
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assumption  that  y  is  not  a  l-CP(i)  with  respect  to  any  constraint  t  6  A.  Note  that  the 
edges  of  the  cone  C\  shown  in  Figure  3(a)  pass  through  the  vertices  of  DUHC[y]  so  that  y 

,  is  almost  -  but  not  quite  -  a  1-ceiling  point  with  respect  to  each  of  the  constraints  binding 

at  x. 

\  To  further  define  the  region  of  only  feasible  integer  solutions,  a  second  “feasibility 

cone”  C2  is  constructed  whose  vertex  is  at  y  and  whose  edges  are  formed  by  constraints 

1 

which  do  not  intersect  UHC[x\.  The  edges  of  C2  that  will  be  used  in  the  proof  and  that  are 
shown  in  Figure  3(a)  represent  a  limiting  case  for  “admissable  non-binding  constraints,” 
i.e.,  constraints  non-binding  at  x  which  neither  chop  off  y  nor  intersect  UHC[x\.  The 
tightest  admissable  non-binding  constraints  must  actually  form  a  cone  strictly  containing 
the  cone  C2  pictured  in  Figure  3(a)  because  the  edges  of  C2  shown  in  the  figure  pass  through 
the  vertices  of  UHC[x ],  i.e.,  they  just  barely  intersect  UHC[x\.  Thus,  integer  points  that 
satisfy  all  of  the  constraints  defining  cones  C\  and  C2  are  feasible  for  all  configurations 
of  the  feasible  region  in  which  y  is  not  a  1-ceiling  point  with  respect  to  any  constraint 
intersecting  UHC[x\. 

Now  a  third  cone  C3  is  defined  such  that  all  points  lying  within  this  cone  have  at 
least  as  good  an  objective  function  value  as  that  of  y  for  all  “admissable  c,”  t.e.,  for  all  c 
such  that  x  is  optimal  for  ( LPr ).  Cone  C3  comes  from  the  Kaxush-Kuhn- Tucker  necessary 
conditions  for  optimality  in  a  linear  program:  at  an  optimal  solution,  the  gradient  of  the 
objective  function  c  can  be  expressed  as  a  nonnegative  combination  of  the  gradients  of 
the  constraints  binding  at  this  solution.  In  other  words,  the  cost  vector  c  lies  in  the  cone 
defined  by  the  row  vectors  of  A  corresponding  to  the  constraints  binding  at  x  [BJ,  p.  215]. 
As  a  result,  an  “optimality  cone”  C3  can  be  constructed  such  that  any  solution  lying  in 
this  cone  possesses  an  objective  value  greater  than  or  equal  to  that  of  y  for  all  admissable 
c.  One  edge  of  C3  is  formed  by  a  constraint  passing  through  y  parallel  to  one  edge  of  C\ 
and  the  other  edge  is  formed  by  another  constraint  passing  through  y  parallel  to  the  other 
edge  of  Ci,  as  illustrated  in  Figure  3(b). 
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Figure  3(a).  General  setting  for  Theorem  4:  Feasibility  Cones  C\  and  C2. 
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Let  region  R  denote  the  set  of  points  that  satisfy  all  six  of  the  constraints  defining 
cones  Ci,C2  C3,  i.e.,  R  =  (C\  D  C2  D  C3).  Since  cones  C2  and  C3  share  a  common 
vertex,  two  of  the  constraints  defining  these  cones  must  be  redundant,  so  that  only  four 
constraints  are  really  needed  to  define  R.  For  example,  in  the  arrangement  shown  in  Figure 
3(b),  cone  C3  appears  to  be  completely  redundant  with  respect  to  C2  (although  this  will 
not  always  be  the  case).  Thus,  all  integer  points  in  R  are  feasible  and  as  good  as  y  for 
all  admissable  c  and  cone  constraints.  We  shall  proceed  to  define  a  region  R  C  R  which 
possesses  only  solutions  strictly  better  than  y. 

First,  we  account  for  the  fact  that  x  may  lie  anywhere  within  UHC[x],  Consider  the 
c<.  ne  C[  generated  by  locating  x  anywhere  strictly  interior  to  UHC[x\  or  on  its  northern  or 
eastern  boundary.  This  cone  is  at  least  as  large  as  another  cone  C"  generated  by  locating 
x  at  some  specific  point  along  the  southern  or  western  boundary  of  U HC\x).  Thus,  it  is 
sufficient  to  consider  the  case  where  x  approaches  some  point  along  the  southern  or  western 
boundary  of  U ffC[x]  because  the  resultant  cone  is  as  small  as  possible.  In  the  following 
analysis,  we  will  examine  the  southern  and  western  boundaries  of  UHC[x\  separately. 

In  Cases  I.  A  and  II.  A,  we  consider  locating  z  simultaneously  at  each  point  along  the 
western  edge  of  U HC\x\,  giving  rise  to  a  sequence  of  cones  as  pictured  in  Figure  3(c).  For 
the  moment,  extend  these  cones  only  as  far  as  vertices  vj  and  v3.  Taking  the  intersection  of 
all  these  truncated  cones  yields  another  truncated  cone  C\  with  extreme  points  v,  tq  and  t>3. 
For  the  case  pictured  in  Figure  3(c),  one  edge  of  cone  C\  passes  through  the  points  u2  and 
V3,  while  the  other  edge  passes  through  the  points  u  1  and  iq.  The  sequence  of  C\  cones,  in 
turn,  gives  rise  to  a  sequence  of  optimality  cones.  In  the  following  figures,  the  cone  labeled 
C3  is  the  intersection  of  all  such  optimality  cones  and  the  region  R  =  (Cj  D  C2  fl  C3 ).  Since 
C2  shares  a  common  vertex  (y)  with  each  optimality  cone  C3,  C2  also  shares  a  common 
vertex  with  C3.  Thus,  only  four  constraints  are  needed  to  define  region  R. 

Recall  that  cTx  >  cTy,  Vx  G  f?,  with  equality  holding  only  when  x  lies  on  the  boundary 
of  cone  C3  and  the  optimal  objective  function  hyperplane  z  =  cTx  coincides  with  an  edge 
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Figure  3(c).  Cone  C\  with  vertex  v  is  the  intersection  of  a  sequence  of  truncated  C\  cones 
whose  vertices  range  between  u\  and  u 2.  Cone  C3  with  vertex  y  is  the  intersection  of  the 
corresponding  sequence  of  C3  optimality  cones. 
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of  C3.  An  edge  of  C3  may  coincide  with  z  =  cTx  only  if  an  edge  of  C\  coincides  with 
z  —  c1  x  because  of  the  way  CVs  edges  are  constructed.  However,  the  corresponding  edges 
of  cones  C\  and  C3  cannot  be  parallel  because  for  every  pair  of  parallel  edges  (one  from 
Ci,  one  from  C3)  at  least  one  of  the  two  edges  becomes  redundant  when  the  intersection 
over  all  Ci  cones  is  taken  to  form  C\  and  the  intersection  over  all  C3  cones  is  taken  to 
form  C3.  Thus,  no  pair  of  parallel  edges  define  R.  Consequently,  cTx  >  cry,\/x  £  R.  even 
if  x  £  R  lies  on  the  boundary  of  the  optimality  cone  C3. 

Cases  I.B  and  II. B  follow  a  plan  similar  to  the  above  except  that  x  is  located  simul¬ 
taneously  at  each  point  along  the  southern  edge  of  UHC\x}.  Taken  together,  these  cases 
will  show  that  a  feasible  integer  solution  x*  which  is  better  than  y  can  be  identified  in  any 
admissible  configuration  of  the  x,  the  binding  constraints,  y  and  c  in  ‘ft2.  To  summarize, 
we  will  contradict  the  optimality  of  y  by  identifying  an  integer  solution  x*  in  each  case 
which  lies  in  the  region  R  defined  by  the  following  four  constraints: 

(In)  and  (lb)  define  the  edges  of  the  truncated  feasibility  cone  C\ , 

(2)  -  defines  one  of  the  two  edges  of  the  feasibility  cone  C2,  and 

(3)  -  defines  one  of  the  two  edges  of  the  optimality  cone  C 3. 

It  is  assumed  that  y  <  x,  but  all  the  arguments  could  be  appropriately  rotated  and  a 
suitable  x*  found,  given  any  other  relationship  between  y  and  .r  ( c.g,,  y  >  x).  Also,  there 
is  no  loss  of  generality  in  letting  y  and  U\  be  defined  as  the  all-integer  vertices  (0.0)  and 
(a,  ?>),  respectively.  Finally,  in  what  follows  the  notation  LHS(j)  will  mean  the  left  hand 
side  of  the  jlh  constraint  evaluated  at  .r*. 
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Case  I. A.  1  <  a  <  b;  x  on  western  edge  of  UHC [x].  See  Figure  4(a). 

( la)  Feasibility  cone  C\  constraint  through  (i>3 ,  u2)  :  (b  +  l)x!  —  (a  —  l)x2  <  b  +  1 

(lb)  Feasibility  cone  C\  constraint  through  (t’i,Uj)  :  —bxt  f  (a  4-  l)x2  <  b 

(2)  Feasibility  cone  C2  constraint  through  (y,  u2)  :  (b  -f  l)xj  —  ax2  >  0 

(3)  Optimality  cone  C3  constraint  through  ( y,  u3 )  :  (b  +  l)xj  —  (a  +  l)x2  <  0 


(i)  a  even:  x*  =  ( I’ll  i  L^J )  =  ( f.  L^J ) 

LHS(la)  =  (b+l)%  -(a-l)L^-J  <  (b  +  1)|  -  (a  -  1)|  =  <  b  <  b  +  1. 

LHS(lb)  —  — 6|  +  (a  +  l)|_^j 


j.  J^+^=§<6,  (beven) 

!  J  |  ^  +  il 


+  >)6  _  b 

=  s±l±i  <  b.  (b  odd) 


LHS(2)  =  (&  4-l)f  >  (6  +  1)|  -  =  0. 

Therefore,  x *  is  feasible. 


LHS( 3)  =  (6  +  l)|-(a  +  l)LS±lJ  <  (6  +  1)|  -  (a  +  1)|  =  ftf*  <  0. 

Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 


(it)  a  odd: 


x*  =  (ffl,f^il)  =  (a±i,r^-1) 


LHS(la)  =  (b+l)l-^-(a-l)\^\  <  (b  +  1)^  -  (a  -  1)^  =  <b+  1. 


LHS(lb)  =  +  (a  +  1)[4±1]  =  | 


=a|=4  +  =  a  +  1  <  6,  (b  even) 

-26-j,  +  {<.+  1)0+1)  =  «±i  <  6.  (b  odd) 


2  1  2 

LHS{2)  =  (6  +  1)^-  -af^-l  >  (6  +  l)l2±H-a^±2i  =  rg±±t»  >  0 

Therefore,  x*  is  feasible. 

TffS(3)  =  (6+l)^-(«  +  l)r^l  =(a+l){S±l-  r^-D^O. 


Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 
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Case  I.B.  1  <  a  <  6;  x  on  southern  edge  of  UHC[x ].  See  Figure  4(b). 

(la)  Feasibility  cone  Cj  constraint  through  («i,U4)  :  —  6i]  +  (a  4-  2)x2  <  6 
(16)  Feasibility  cone  C\  constraint  through  (u3,ui) :  6xj  -  (a  -  l)x2  <  b 
(21  Feasibility  cone  C2  constraint  through  (y,u4)  :  — 6xi  +  (a  +  1)^2  >  0 
(3)  Optimality  cone  C3  constraint  through  (y,U])  :  bxx  -  ax2  >  0 

(i)  a  even,  a  <  6  :  x*  =  (|,  [|j) 

LHS(la)  =  -6(f)  +  (a  +  2)[|j  <  -6(f)  +  (a  +  2)f  =  6. 

LHS(lb)  =  6(f)  -  (a  -  1)L|J  <  f  -  («  -  <  b- 

LHS(2)  =  -6(f)  +  (a  +  l)[fj  >  -f  +  (a  +  l)*f±  =  >  o. 

Therefore,  x*  is  feasible. 

I/f5(3)  =  6(f)  — aLfj  >a(|-  lfj)>0. 

Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 


(ii)  a  odd, a  <  6  : 

LHS(la)  —  — 6^,a^1^  +  (a  +  2)f 


x*  =  ( 


II  < 


—  ft(a  + 1 )  . 

2 

-6(a  +  l)  . 
2  ‘r 


(a+2)(6+l)  _  a+b+2  <  l 
2  '2  —  1 
(o+2)6  _  b  <  l 
2  —  2  — 


(b  odd) 
(b  even) 


LHS(lb)  =  b~^-  -  (a  -  l)f|l  <  ^  “  (a  ~  l)f  =  6. 


LHS{2)  =  -6^  +  (a  +  l)[fl  =  (a  +  l)(f£]  -  f )  >  0. 


Therefore,  x*  is  feasible. 

LHS( 3)  =  6^  -  affl  >  (a  +  l)f  -  a^)  =  ^  >  0. 
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Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 

(■“)»  = *'  =  (lfj.l£j) 

IHS(U)  =  -Hf J  +  (a  +  2)lij  =  -H|J  +  (b  +  2)[jJ  =  2L|j  <  b. 

LHS(lb)  =  HfJ  -  (a  -  1)L}J  =  [jJ  <  t. 

lh s(2)  =  -HU  +  (a  +  1)L|J  =  -HtJ  +  (1  +  i)L|J  =  LlJ  >  o. 

Therefore,  x*  is  feasible. 

tHS( 3)  =  HU  -  «lU  »  0. 

Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 


4.  Search  Constraints 


Case  II.A.  1  <  6  <  a;  x  on  western  edge  of  UHC[x\.  See  Figure  4(c). 

(la)  Feasibility  cone  C\  constraint  through  (i>4,  u2) :  (6  +  2)xj  —  ax2  <  a 
(16)  Feasibility  cone  C\  constraint  through  (i>2,tii) :  —('■>  —  l)xj  +  ax2  <  a 

(2)  Feasibility  cone  C2  constraint  through  (y,  u2)  :  (6  +  l)xi  —  ax2  >  0 

(3)  Optimality  cone  C3  constraint  through  (y,ui) :  6xj  —  ax2  <  0 

(i)  6 even  :  x*  =  (IfJ,  1§J)  = 

LHS(la)  =  (6  +  2)L|j  -  f  <  a^)  -  f  =  a. 

LHS(lb)  =  -(6  -  l)LfJ  +  «(|)  <  -(6  -  l)2fi  +  f  f  =  a. 

Ltf5(2)  =  (6  +  l)LfJ  ~^>(6  +  l)f +  f  =  ^<0. 

Therefore,  x*  is  feasible. 

LffS(3)  =  6Lf J  -  f  =  KLfJ  -  f)  <  0. 

Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 

(ii)  ft  odd  :  i*  =  ( Tfl.liJ)  =  (ffl.  |) 

LHS(\a)  =  (t  +  2)f|l  -  =  (6  +  l)(f|l  -  |)  +  ffl  =  {  Lfu  <  £  ^d” 

LHS(lb)  =  —(6  -  l)r|l  -  o(|)  <  -J(6-  1)  +  $  <  «■ 

LHS{2)  =  (6  +  l)r|l  -  >  (i  +  1)|  -  =  0. 

Therefore,  x*  is  feasible. 

LHS( 3)  =  6f|l  -  <  6^  “  =  =$**  <  0. 

Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 
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Case  II.B.  1  <  6  <  a;  i  on  southern  edge  of  UHC[x].  See  Figure  4(d). 

(la)  Feasibility  cone  C\  constraint  through  (t>2»“4) :  —(6  —  l)*i  +  (a  +  l)x2  <  a  +  1 
(16)  Feasibility  cone  C\  constraint  through  (t>4,ui) :  (6  +  l)xj  —  ax%  <  a 

(2)  Feasibility  cone  C2  constraint  through  (y,  u4)  :  -bxx  +  (a  +  l)x2  >  0 

(3)  Optimality  cone  C3  constraint  through  (y, u3) :  (6  H  l)xj  -  (a  +  l)x2  >  0 

(i)  6 even  :  x*  =  (L^J,  Til)  =  f) 

LHS(\a)  =  -(6  -  l)[2flj  +  <  -(6  -  1)|  +  <  a  +  1. 

X^S(16)  =  (6+l)[a±iJ  -a(|)  <  (a  +  l)^±i-  f  =  a±|±i  <  a. 

LHS( 2)  =  >  -6^  +  =  o. 

Therefore,  x*  is  feasible. 

LHS( 3)  =  (6  +  l)ia±lj  -  >  (6+  l)f  -  >  0. 

Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 

(n)  6odd :  rfi)  = 

Iff5(la)  =  -(6  -  l)fsf*-l  +  ^^±12  <  -(6  -  l)i^  +  fr+W+l)  =  a  +  1. 
IJf5(16)  =  (6+  l)f*fi)  -  2^  <(6+l)i2±21-(^)  =  6  +  l  <a. 

LHS( 2)  =  -6fs±2-1  +  >  -6&±21  +  ls±lK*±U  =  «=i±I  >  0. 

Therefore,  x*  is  feasible. 

LHS{ 3)  =  (6  +  ljfs^-l  -  ^±1K*±H  >  (fc  +  i)ii±H  -  =  0. 

Thus,  x*  is  strictly  better  than  y,  contradicting  the  assumption  that  y  is  optimal. 
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The  final  case  that  needs  to  be  considered  is  where  either  a  =  0  or  b  =  0  (but  not  both). 
If  a  =  0,  the  integer  solution  x*  =  v2  =  (0, 1)  is  clearly  feasible  and  better  than  y,  while  if 
6  =  0,  the  same  can  be  said  of  x*  ~  v3  —  (1,0).  I 

Conjecture  4a  (Theorem  4  with  n  >  3)  Suppose  an  ( ILP )  is  non-empty  and  bounded. 
Further  suppose  that  c  ^  0  and  x  is  an  optimal  solution  for  (LPr).  Then  every  optimal 
solution  for  (ILP)  is  a  feasible  1-CP(*‘),  where  constraint  (i)  intersects  each  of  the  2“ 
UHC[x\ s. 

Rationale:  The  reasoning  behind  this  is  that  it  seems  the  “geometry”  of  the  situation 
in  higher  dimensions  will  reflect  that  of  the  two-dimensional  case  described  in  the  proof  of 
Theorem  4.  For  instance,  projecting  the  situation  in  R3  onto  any  of  the  two-dimensional 
coordinate  planes  xi-x2,  x2-x3,  or  Xj-x3  would  reveal  a  feasible  integer  solution  x*  in  this 
projection  with  a  value  superior  to  that  of  y.  Let  p  denote  the  midway  point  between 
x  and  y  and  U HC\p ]  the  unit  hypercube  containing  p  with  all-integer  vertices.  For  each 
case  of  Theorem  4,  the  solution  x*  is  the  feasible  integer  solution  closest  to  p,  which  is 
where  the  region  defined  by  the  three  cones  is  widest.  In  higher  dimensions,  it  seems  likely 
that  among  the  lattice  points  between  y  and  UHC[x\,  one  such  integer  solution  which 
is  approximately  halfway  between  x  and  y  will  satisfy  all  the  requirements  for  a  feasible 
solution  with  an  objective  value  greater  than  y  for  all  admissable  c,  i.e.,  all  c  such  that 
x  is  optimal  for  ( LPr ).  First  note  that  for  some  vertex  v  of  UHC\p],  every  component  vj 
of  this  vertex  differs  from  p;  by  no  more  than  i.  Then  observe  that,  by  construction,  the 
width  of  the  cone  Cj  when  sliced  through  p  parallel  to  each  coordinate  axis  is  at  least  1. 
Furthermore,  the  width  of  both  cones  C2  and  C3  when  sliced  through  p  parallel  to  each 
coordinate  axis  is  at  least  j.  So,  in  the  neighborhood  of  p,  it  appears  possible  to  move 
to  an  integer  solution  x*  (a  vertex  of  U HC\p])  which  satisfies  the  constraints  defining  all 
three  cones.  Thus,  it  seems  plausible  in  higher  dimensions  that  we  will  be  able  to  find  a 
feasible  integer  solution  x*  better  than  y  when  we  assume  that  y  is  not  a  l-CP(i)  with 
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respect  to  any  constraint  which  intersects  one  the  2“  £//fC[x]’s.  I 

Consider  the  implications  for  developing  an  algorithm  to  find  an  optimal  solution 
for  ( ILP )  if  this  conjecture  is  correct.  The  conjecture  implies  that  the  set  of  functional 
constraints  whose  l-CP(»)’s  we  must  enumerate  in  order  to  solve  ( ILP )  consists  of  just 
those  which  intersect  a  unit  hypercube  with  all- integer  vertices  containing  x.  Part  of  this 
set  of  constraints,  namely,  the  set  A  of  constraints  binding  at  x,  may  be  identified  by  using 
the  simplex  method  to  solve  the  ( LPr ).  The  next  step  would  be  to  pivot  from  x  to  find 
all  extreme  points  adjacent  to  it  and  check  whether  or  not  each  adjacent  extreme  point  is 
contained  in  UHC[x\,  if  it  is,  then  any  constraint  binding  at  this  extreme  point  intersects 
UHC[x].  We  still  would  need  to  test  whether  or  not  any  of  the  remaining  constraints 
intersect  UHC[x].  A  sufficient  condition  for  excluding  a  constraint  from  the  set  is  if  the 
distance  from  the  constraint  to  x  at  its  closest  point  is  more  than  \/n.  Such  a  constraint 
cannot  possibly  intersect  UHC[x\.  Constraints  less  than  a  distance  of  y/n  from  x  at 
their  closest  point  do  not  neccessarily  intersect  UHC[x],  but  to  determine  this  precisely 
takes  more  effort.  However,  computational  experience  has  shown  that  it  is  unnecessary 
to  determine  the  complete  set  of  constraints  intersecting  UHC[x\  because  the  problem 
is  usually  solved  by  the  time  we  have  enumerated  the  1-ceiling  points  with  respect  to 
just  a  few  constraints  binding  at  x.  Thus,  while  Theorem  4  and  Conjecture  4a  are  of 
theoretical  interest,  they  have  had  little  impact  on  the  development  of  the  Exact  Ceiling 
Point  Algorithm  described  in  a  subsequent  report. 


5.  Underlying  Assumptions 


The  theorems  of  the  previous  sections  imply  that  to  solve  a  non-trivial,  non-empty, 
bounded  ( ILP )  with  c  ^  0,  it  is  sufficient  to  enumerate  its  feasible  l-CP(*)’s.  In  this 
section  we  consider  the  effect  of  placing  these  conditions  on  (ILP).  First,  assuming  that 
c  ^  0  just  means  that  we  have  an  optimization  problem.  Second,  to  check  that  (ILP)  is 
not  trivial  to  solve,  we  simply  examine  some  of  the  data;  if  both  c  <  0  and  b  >  0,  then  i  =  0 
trivially  solves  the  (ILP).  Third,  although  we  can  conclude  that  (ILP)  is  empty  if  the 
associated  (LPr)  is  empty,  the  contrapositive  is  not  true.  So,  when  ( LPr )  is  non-empty, 
we  can  determine  whether  or  not  (ILP)  is  empty  by  seeking  only  l-CP(FR)'s  because  an 
(ILP)  containing  no  1-ceiling  points  does  not  contain  any  feasible  integer  solutions  either, 
and  vice-versa,  as  Theorem  5  demonstrates.  Theorem  5  follows  directly  from  the  next  two 
lemmas. 

Lemma  5a  The  set  of  feasible  integer  solutions  for  (ILP)  is  empty  if  and  only  if  the  set 
of  CP(FR)' s  is  empty. 

Proof:  (=>)  Assume  no  CP(FR)' s  exist.  Then,  by  Lemma  la,  no  extreme  points  of 
the  convex  hull  of  feasible  integer  solutions  exist.  Since  the  feasible  region  for  (ILP)  is 
bounded  below,  the  convex  hull  of  feasible  integer  solutions  must  be  empty  (rather  than 
being  unbounded).  Therefore,  no  feasible  integer  solutions  exist. 

(•$=)  Assume  no  feasible  integer  solutions  exist.  Then,  since  every  CP(FR)  is  a  special 
kind  of  feasible  integer  solution,  no  CP(FR)’ s  may  exist  either.  I 

Lemma  5b  The  set  of  CP(FR)’ s  for  (ILP)  is  empty  if  and  only  if  the  set  of  \-CP(FR)'s 
is  empty. 

Proof:  (=*■)  Assume  no  l-CP(Ffl)’s  exist.  Since  every  CP(FR)  is  a  special  kind  of 
feasible  1  -CP(FR)  by  Lemma  lc,  part  (a),  no  CP(FF)’s  may  exist  either. 
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(-4=)  Assume  no  CP(FR)'s  exist.  Then,  by  Lemma  5a,  the  set  of  feasible  integer  solutions 
for  ( ILP )  is  empty  and  thus  no  \-CP(FR)'s  exist.  I 

Theorem  5  The  set  of  feasible  integer  solutions  for  (ILP)  is  empty  if  and  only  if  the  set 
of  l-CP(FR)'s  is  empty. 

Proof:  Immediate  from  Lemmas  5a  and  5b.  I 

Finally,  the  boundedness  condition  can  also  be  checked  by  solving  the  linear  program¬ 
ming  relaxation  of  (ILP)  and  possibly  enumerating  1-ceiling  points. 

Theorem  6  If  the  objective  function  of  (LPr)  is  bounded  above,  then  either  the  objective 
function  of  (ILP)  is  also  bounded  above  or  (ILP)  is  infeasible. 

Proof:  Even  if  solutions  exist  for  (LPr),  infeasiblity  of  (ILP)  is  clearly  a  possibility. 
For  example,  the  region  defined  by  {.r|  1  <  3.r;  <  2 ,Vj}  contains  real  but  no  integer 
solutions.  Because  (LPr)  is  a  relaxation  of  (ILP)  the  optimal  objective  function  value  of 
(ILP)  must  be  less  than  or  equal  to  that  of  (LPr).  Thus,  if  the  former  is  bounded  above, 
the  latter  must  also  be  bounded  above.  I 

When  the  objective  function  of  (LPr)  is  unbounded  and  all  data  are  rational,  we  can 
use  the  following  result  (see  [Pa,  Lemma  3]  or  [BG,  Theorem  1]). 

Theorem  7  (Papadimitriou)  Assume  the  data  of  (ILP)  are  rational.  If  the  objective 
function  of  (LPr)  is  unbounded  above,  then  either  the  objective  function  of  (ILP)  is  also 
unbounded  above  or  (ILP)  is  infeasible. 

Whether  (LPr)  is  bounded  or  not,  we  may  use  an  algorithm  which  solves  the  linear 
programming  relaxation  and  searches  only  for  1-ceiling  points  of  (ILP)  to  make  a  relevant 
statement  about  the  optimal  value  of  (ILP).  To  summarize:  (1)  if  (LPr)  is  infeasible. 
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so  is  ( ILP );  (2)  if  the  objective  function  of  ( LPr )  is  bounded  above,  then  a  1-ceiling 
point  algorithm  applied  to  (ILP)  will  determine  either  that  (ILP)  has  no  feasible  integer 
solutions  or  that  it  also  has  an  objective  function  which  is  bounded  above;  (3)  if  the 
objective  function  of  (LPr)  is  unbounded,  then  a  1-ceiling  point  algorithm  applied  to 
(ILP)  will  conclude  either  that  (ILP)  has  no  feasible  solutions  or,  upon  first  identifying 
a  feasible  integer  solution,  that  (ILP)'s  objective  function  is  also  unbounded  above.  In 
this  latter  case,  the  integer  linear  program  has  most  likely  been  formulated  incorrectly. 

6.  Summary  of  Concepts  with  an  Example 

We  shall  now  briefly  review  the  ideas  presented  in  this  report  and  try  to  illustrate  their 
effects  graphically.  The  key  theorems  given  above  will  be  applied  in  turn  to  the  simple 
example  of  Figure  1  having  two  variables  and  three  constraints.  After  each  is  applied,  the 
set  of  integer  solutions  on  which  we  need  to  focus  in  order  to  solve  the  problem,  shown  as 
the  solid  lattice  points  in  Figures  5  and  6,  is  reduced.  The  first  diagram  below,  Figure  5(a), 
illustrates  that  initially  all  feasible  integer  solutions  of  our  (ILP)  are  under  consideration. 
The  significance  of  Theorem  2  is  that  to  solve  our  problem,  we  need  only  to  focus  on  the 
feasible  l-CP(i)’s,  as  shown  in  Figure  5(b). 

Considering  Theorem  3,  we  need  retain  from  the  set  of  all  feasible  l-CP(t),s  only  those 
which  are  1-ceiling  points  with  respect  to  afunctional  constraint.  Therefore,  we  can  ignore 
those  integer  points  which  lie  along  either  of  the  axes  but  which  are  either  infeasible  or 
are  not  1-ceiling  points  with  respect  to  a  functional  constraint.  The  effect  of  this  is  shown 
in  Figure  6(a),  where  the  set  of  1-ceiling  points  {(xi,0),Xi  =  0, 1,...,4}  are  dropped  from 
consideration. 

Finally,  if  we  solve  the  linear  programming  relaxation  associated  with  (ILP),  we  may 
take  one  more  step  based  on  its  optimal  solution,  x.  From  the  remaining  set  of  points, 
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i  i 

Figure  5. (a)  All  feasible  solutions.  (b)  Feasible  l-CF(i)’s  (Theorem  2). 

Theorem  4  allows  us  to  just  concentrate  on  those  which  are  1-ceiling  points  with  respect  to 
a  functional  constraint  that  intersects  a  unit  cube  containing  x  and  possessing  all-integer 
vertices.  If  we  assume  that  x  occurs  at  the  intersection  of  constraints  (1)  and  (2),  then  we 
no  longer  need  to  examine  most  of  the  feasible  l-CP(3)’s,  namely,  {(4,  x2 ),  x2  =  0, 1,2}, 
because  they  are  not  l-CP(i)’s  with  respect  to  either  constraint  (1)  or  constraint  (2).  The 
remaining  integer  points  which  are  still  sufficient  to  solve  ( ILP )  are  shown  in  Figure  6(b). 

Thus,  the  net  effect  of  the  theorems  is  to  significantly  reduce  the  set  of  feasible  integer 
solutions  on  which  we  need  to  focus  our  attention  in  order  to  solve  (ILP),  assuming  that 
feasible  solutions  exist.  Subsequent  reports  will  describe  a  heuristic  algorithm  and  an  exact 
algorithm,  respectively,  which  are  designed  to  take  advantage  of  these  ideas. 
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This  report  examines  the  role  played  by  several  kinds  of  ceiling  points  in  solving  the  pure, 
general  integer  linear  programming  problem  (ILP).  While  no  assumptions  are  made 
concerning  the  structure  or  signs  of  the  data  of  the  problem,  it  is  assumed  that  the  feasible 
region  for  (ILP)  is  non-empty  and  bounded.  A  ceiling  point  with  respect  to  a  single 
constraint  may  be  thought  of  as  an  integer  solution  on  or  close  to  the  boundary  of  the 
feasible  region  defined  by  the  constraint.  The  definition  of  a  ceiling  point  with  respect  to  a 
single  constraint  is  extended  to  take  multiple  constraints  into  consideration  simultaneously, 
defining  what  is  called  a  feasible  ceiling  point  It  is  shown  that  the  set  of  all  feasible  ceiling 
points  contains  at  least  one  optimal  solution  for  (ILP).  A  related  class  of  solutions  called 
feasible  1-ceiling  points  is  also  characterized  and  shown  to  contain  all  optimal  solutions  for 
(ILP).  Moreover,  1 -ceiling  points  are  computationally  easier  to  identify  than  ordinary  ceiling 
points  and  may  be  sought  with  respect  to  one  cons  taint  at  a  time.  It  is  also  demonstrated 
that  solving  (ILP)  requires  only  enumerating  feasible  1 -ceiling  points  with  respect  to  a 
subset  of  all  functional  constraints. 
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